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We study a systematic improvement of perturbation theory for gauge fields on the lattice [ffl]; the 
improvement entails resumming, to all orders in the coupling constant, a dominant subclass of 
tadpole diagrams. 

This method, originally proposed for the Wilson gluon action is extended here to encompass 
all possible gluon actions made of closed Wilson loops; any fermion action can be employed as 
well. The effect of resummation is to replace various parameters in the action (coupling constant, 
Symanzik and clover coefficient) by "dressed" values; the latter are solutions to certain coupled 
integral equations, which are easy to solve numerically. 

Some positive features of this method are: a) It is gauge invariant, b) it can be systematically 
applied to improve (to all orders) results obtained at any given order in perturbation theory, c) it 
does indeed absorb in the dressed parameters the bulk of tadpole contributions. 
Two different applications are presented: The additive renormalization of fermion masses, and 
the multiplicative renormalization Zy (Z^) of the vector (axial) current. In many cases where 
non-perturbative estimates of renormalization functions are also available for comparison, the 
agreement with improved perturbative results is consistently better as compared to results from 
bare perturbation theory. 
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1. Introduction 



It is well known that quantities measured through numerical simulation are characterized by 
significant renormalization effects, which must be properly taken into account before making any 
comparisons to corresponding physical observables. 

Although the renormalization procedure can be formally carried out in a systematic way to any 
given order in perturbation theory, calculations are notoriously difficult, as compared to continuum 
regularization schemes. Furthermore, the convergence rate of the resulting asymptotic series is 
often unsatisfactory. 

Some years ago, a method was proposed to sum up a whole subclass of tadpole diagrams, 
dubbed "cactus" diagrams, to all orders in perturbation theory [Q, This procedure is gauge 
invariant, it can be systematically applied to improve (to all orders) results obtained at any given 
order in perturbation theory, and it does indeed absorb the bulk of tadpole contributions into an in- 
tricate redefinition of the coupling constant. The agreement of available non-perturbative estimates 
of renormalization coefficients with cactus improved perturbative results is consistently better as 
compared to results from bare perturbation theory. 

In the present work we extend the improved perturbation theory method of Refs. [|2|, Q], to 
encompass the large class of actions (including Symanzik improved gluon actions combined with 
any fermionic action) which are used nowadays in simulations of QCD. In Section II we present 
our calculation, deriving expressions for a dressed gluon propagator. The methodology can be also 
applied to dress the gluon and fermion vertices (appears in Ref. [[]]]). We show how these dressed 
constituents are employed to improve 1-loop and 2-loop Feynman diagrams coming from bare 
perturbation theory. In Section III we apply our improved renormalizaton procedure to a number 
of test cases involving Symanzik gluons and Wilson/clover/overlap fermions. 



2. The Method 



In this Section, we start illustrating our method by showing how the gluon propagator is 
dressed by the inclusion of cactus diagrams. We will then explain how this procedure is applied to 
Feynman diagrams at a given order in bare perturbation theory, concentrating on the 1- and 2-loop 
case. 

2.1 Dressing the propagator 

We consider the Symanzik improved gluon action involving Wilson loops with up to 6 links: 



Sg — —j 
So 



C Re Tr ( 1 ~ ^plaquette ) + C\ £ Re Tr ( 1 - Rectangle ) 

plaquette 



rectangle 



C 2 ^ Re Tr ( 1 - U chail ) + C 3 £ Re Tr ( 1 - ^/parallelogram ) 
chair parallelogram 



(2.1) 



The coefficients a satisfy a normalization condition co + 8ci + I6C2 + 803 = 1 which ensures the 
correct classical continuum limit of the action. 
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The quantities [/, (i = O(plaquette), 1 (rectangle), 2 (chair), 3 (parallelogram)) in Eq. ( |2T| ) are 
products of link variables U Xill around the perimeter of the closed loop. Using the Baker-Campbell- 
Hausdorff (BCH) formula, Uj takes the form: 

Ui = exp (igoF^ + iglF® + iglF® + tffejj)) (2.2) 

where F>^ is simply the sum of the gauge fields on the links of loop i, while Ff^ (j > 1) are j-th 
degree polynomials in the gauge fields, constructed from nested commutators. 



Let us define the cactus diagrams which dress the gluon propa- 
gator: These are tadpole diagrams which become disconnected 
if any one of their vertices is removed; further, each vertex 
is constructed solely from the Zy parts of the action. A di- 
agrammatic equation for the dressed gluon propagator (thick 
line) in terms of the bare propagator (thin line) and 1 -particle 
irreducible (1PI) vertices (solid circle) reads: 




Fig.l: A cactus 



The 1PI vertex obeys the following recursive equation: 



(2.3) 



+ do 



(2.4) 



The bare inverse gluon propagator D 1 results from the total gluon action ([[lj]) and reads 

D-^(k)=Y t (P p S^-k^k p S pv )d» p + ^= £ ( CiG %(k)\+hbL (2.5) 

O K ' 1 S i=0.1.2.3 V / 1 b 



where Co = Co + 8c i + 16c2 + 8C3, C\ = C2 + C3, C2 = c\ — C2 — C3. For further definitions 
of quantities appearing above, the reader can refer to [[]]]. The matrices G^(k) are symmetric and 
transverse, and originate from a Tr (/y Zy ) term of the gluon action. Consequently, the diagrams 
on the r.h.s. of Eq. (^4), are a linear combination of G^ l >(k); this implies that the 1PI vertex G m (k) 



(the l.h.s. of Eq. (|2.4|)) can be written as: 



G w \k) = Ob G (0) (Jfc) + «! (it) + a 2 G (2) (k) + a 3 G (3) (k) 



(2.6) 



Each of the quantities a, will in general depend on N, go, Co, ci , C2, C3, but not on the momentum. 
Eq. (2.3) leads to the following expression for the inverse dressed propagator (D dr )~ 1 (k) 

1 



(Z)*) " 1 = co G [U > + cj G l ^ + c 2 + c 3 ^ + 



:(0) 



(i). 



;(2) 



:(3). 



kflky . 



Ci = Ci 



a, 



(2.7) 
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We observe that dressing replaces the bare coefficients c, with improved ones c,- , and leaves 
the longitudinal part intact. This property ensures gauge invariance of the results. 
In terms of the dressed propagator, Eq. ( |2.4[ ) can be drawn as: 



(2.8) 



A typical diagram on the r.h.s. of Eq. (2JS) is the sum of 4 terms, and has (n — 2) /2 1-loop integrals 
in the diagram (coming from the contraction of two powers of via a dressed propagator), and 
will contribute one power of j8,-(co,c 1,02,03), where: 

r* d 4 q 



ft) 
ft 
ft 
ft 



2 ql D* (<?) -2q^q v £>J V (<?)) 

(4^-§*)D* (?) + ^(4-^)D*(?)-2? M ? v (4-^)Df v (<?)) 



-» (2tt) 4 
-* (2tt) 4 

1 (2^4 ( 8 -^vXp (9)/2 " ^ <?P (8-? 2 v)<p (<?)/2 



t p7r)J v 3^(4-^ v )£>J p ( ? )/2-3^^(4-^)Dj v ( ? )/2 



(2.9) 



(jU,v,p assume distinct values; no summation implied). We note that j3; are gauge independent, 
since the longitudinal part cancels in the loop contraction. 



In order to set Eq. (2.4) in a mathematical form, we need to evaluate F(n;N) which is the 
sum over all complete pairwise contractions of Tr{T ai T" 2 . . . T a " }. Use of F(n;N), along with the 
integrals ( |2.9[ ), allows us to resum (2.4), leading to 



(N 2 -l) 



-p,gl(N-l)/(4N) 



N-l 
N 



4-i(*o ft/2) +24- 2 (So ft/2) (2-10) 



(Ljj(x): Laguerre polynomials). Eqs. ( 2.10 ) are 4 separate equations where unknown quantities are 
the coefficients a,- ; they appear on the l.h.s., as well as inside the integrals j8; of the r.h.s, by virtue 
of Eqs. (2^, 2.7). It is worth mentioning that all combinatorial weights are correctly incorporated 
in the procedure. 

Eqs. dUOb can be solved numerically and each choice of values for (c,- , go , AO leads to a set of 
values for c, = c, — a, that are no longer normalized; one may express the results of our procedure 
in terms of a normalized set of improved coefficients, c,/Cb and an improved coupling constant 
= gl/Co, where: Co = cq + 8ci + I6C2 + 8C3 . For reasons of simplicity we define rescaled 
quantities: 



Ci 

Yi = — 
80 



8l 



ft(co,ci,c 2 ,c 3 ) =goPi(co,ci,c 2 ,c 3 ) = Pi(%,fi,%Y}) (2.11) 



fi must now satisfy the coupled equations: 
1 



N 2 -l 



Tie 



A ^- 1)/(W) '— 4-i(A/2)+24-2(A/2) 



N 



(2.12) 



For SU (2) and SU (3) the Laguerre polynomials have a simple form and Eqs. (2.12) can be written 
explicitly ||Xjj . 
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Since Eqs. (2.12) have the form x = f(x), they can be numerically solved using a fixed point 
procedure [jlj|. A unique solution for fj always exists for all physically interesting values of c,, and 
for all values of go well inside the strong coupling region. The convergence of the procedure has 
been verified in a number of extreme cases. 



2.2 Numerical values of improved coefficients 

We now present the values of the dressed coefficients for several gluon actions of interest. In 
Figs. 2-5, one can see the improved coefficients for Plaquette, Tree-level Symanzik, Iwasaki and 
Tadpole improved Liischer-Weisz actions. Results for DBW2 action are listed in Table I. 



UNO 
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Fig. 2: Improved coefficient Co for N—2 and A^=3 
(plaquette action) 
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Fig. 3: Improved coefficients cq and ci (tree-level 
Symanzik improved action, N = 3) 
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Fig.4: Improved coefficients Co and ci (Iwasaki 
action, N = 3) 



Fig. 5: Coefficients c,- and their dressed counterparts c ( - for 
different values of j3 cq = 6co/^q (TILW actions, N = 3) 



P = 6/ 'si 


CO 


c\ 


C() 


Cl 


1.1636 


5.29078 


-0.53635 


3.39826 


-0.22528 


0.6508 


12.2688 


-1.4086 


8.8070 


-0.7313 



TABLE I. Improved coefficients Co and 5i in the DBW2 action (cq and c\ are 
obtained starting from j3 cq = 6.0 and 6.3) 



3. Applications 

We now turn to two different applications of cactus improvement: The additive mass renor- 



5 



Cactus Improvement 



Martha Constantinou 



malization for clover fermions and the 1-loop renormalization of the axial and vector currents using 
the overlap action. Both cases employ Symanzik improved gluons; hence, our results are presented 
for various sets of Symanzik coefficients. 

3.1 Critical mass of clover fermions 

It is well known that an ultra-local discretization of the fermion action without doubling breaks 
chirality. Consequently, we must demand a zero renormalized fermion mass, in order to ensure chi- 
ral symmetry while approaching the continuum limit. For this purpose, the bare mass is additively 
renormalized from its zero tree-level value to a critical value dm. 

We calculate the 1-loop result for the critical mass dm\-\ oop using clover fermions and Symanzik 
improved gluons. The result is then dressed with cactus diagrams in order to get the improved value 
d ,n i-ioop- Details on the definition of dm as well as a 2-loop calculation of dm with the same ac- 
tions can be found in Ref. [|6],f7],[8|]. The result of the 1-loop diagrams contributing to <imi_i 00 p can 
be written as a polynomial in the clover parameter, and is independent of the number of fermion 
flavors Nf. Some numerical values for <iwii_i oop corresponding to the plaquette and Iwasaki actions 
(N = 3) appear in Ref. [jl|]. 

Using the critical mass, one can evaluate the critical hopping parameter, K^ r = 1/ (2dm + 8r) 
(r is the Wilson parameter). Estimates of fQ r from numerical simulations exist in the literature 
for the plaquette action |§j(A^ = 0), [JTTJ] (Nf = 2), and also the Iwasaki action ||T|] (Nf = 2). 



Perturbative (unimproved and dressed) and non-perturbative results are listed in Table II for specific 
values of csw- It is clear that cactus dressing leads to results for k cy which are much closer to values 
obtained from simulations. 



Action 


N f 




csw 


cr,l— loop 


cr. 1— loop 


non— pert 
Kcr 


Plaquette 





6.00 


1.479 


0.1301 


0.1362 


0.1392 


Plaquette 





6.00 


1.769 


0.1275 


0.1337 


0.1353 


Plaquette 


2 


5.29 


1.9192 


0.1262 


0.1353 


0.1373 


Iwasaki 


2 


1.95 


1.53 


0.1292 


0.1388 


0.1421 



TABLE II. 1-loop results and non-perturbative values for K a 

3.2 One-loop renormalization of fermionic currents 

As a second application of cactus improvement, we investigate the renormalization constant Zy 
(Z A ) of the flavor non-singlet vector (axial) current in 1-loop perturbation theory. Overlap fermions 
and Symanzik improved gluons are employed. Bare 1-loop results for Zva have been computed in 
the literature [ jT3| , §, they depend on the overlap parameter p (0 < p < 2). 



The renormalization constants Zy and Z A are equal [ |13| ] when using the overlap action and 
in the MS scheme. Table III of Ref. [jl|] gives the values of Zva and Zy A for different sets of the 
Symanzik coefficients, choosing p = 1.0, p = 1.4. The dependence of Zy^ and Zy A on the overlap 
parameter p is shown in Fig. 6, where we plot our results for three actions: Plaquette, Iwasaki and 
TILW. Note that improvement is more apparent for the case of the plaquette action. This was to 
be expected, since improved gluon actions were constructed in a way as to reduce lattice artifacts, 
in the first place. A comparison between our improved Zy^ values and some non-perturbative 
estimates jjTBj], shows that improvement moves in the right direction. 
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Fig. 6: Plots of Zy& and Zy A for the plaquette, Iwasaki and TILW actions. Labels have been placed in the 
same top-to-bottom order as their corresponding curves. 
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